Ab Initio Molecular Dynamics Study of Aqueous Solvation of Ethanol and Ethylene 
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The structure and dynamics of aqueous solvation of ethanol and ethylene are studied by DFT- 
based Car-Parrinello molecular dynamics. We did not find an enhancement of the structure of 
the hydrogen bonded network of hydrating water molecules. Both ethanol and ethylene can easily 
be accommodated in the hydrogen-bonded network of water molecules without altering its struc- 
ture. This is supports the conclusion from recent neutron diffraction experiments that there is no 
hydrophobic hydration around small hydrophobic groups. Analysis of the electronic charge dis- 
tribution using Wannier functions shows that the dipole moment of ethanol increases from 1.8 D 
to 3.1 D upon solvation, while the apolar ethylene molecule attains an average dipole moment of 
0.5 D. For ethylene, we identified configurations with tt-H bonded water molecules, that have rare 
four-fold hydrogen-bonded water coordination, yielding instantaneous dipole moments of ethylene 
of up to 1 D. The results provide valuable information for the improvement of empirical force fields, 
and point out that for an accurate description of the aqueous solvation of ethanol, and even of the 
apolar ethylene, polarizable force fields are required. 
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I. INTRODUCTION 

The study of the solvation of alcohols in aqueous so- 
lution is of fundamental interest in physics, biology and 
chemistry, but also of importance in technical applica- 
tions § 1 1 1, 1 1 1, 1 1 0. Among the alcohols, 
ethanol is one the most well studied compounds. Aque- 
ous ethanol solutions are a common in chemical research 
and industry applications. Ethanol can be produced by 
aqueous hydration of ethylene, that is readily available 
from natural sources. This process that can be acceler- 
ated by acid catalysis. The reverse route of decomposing 
of ethanol into water and ethylene is also of importance. 
For example, for developing countries who do not have 
a large supply of fossil fuels, dehydration of ethanol ob- 
tained from biomass is often the most economical way 
to produce ethylene. As is well known, solvation struc- 
tures play a crucial role in aqueous solution chemistry, 
where reactive events often require a significant reorder- 
ing of the water molecules in the solvation shell. The 
solvation of ethanol and ethylene are therefore crucial in 
the course of their (acid-catalyzed) interconversion. The 
aqueous solvation of these molecules will be addressed in 
the present paper. Elsewhere we report on an ab initio 
molecular dynamics study of the acid-catalyzed conver- 
sion 

Ethanol and ethylene have distinct solvation proper- 
ties in aqueous solution. Ethanol is easily soluble as 
its polar hydroxyl group can participate in the hydrogen 
bonded network and the hydrophobic ethyl group is rel- 
atively small. In contrast, the apolar ethylene molecule 
has a much weaker interaction with water and is gener- 
ally considered to be hydrophobic. Mixtures of water and 
ethanol have been studied extensively, both experimen- 
tally as by molecular simulation. Experimental studies 
employing NMR (l3[ ^ , ultrasonic absorption 

p8[ , infrared absorption spectroscopy |l5|, |l9|, pO[ , 
mass spectroscopy [G^ Ell, X-ray diffraction measure- 




and di- 
81 have 



ments neutron diffraction ||22| , [2E 

electric relaxation measurements |25|7p6[ 
been performed to unravel the solvation properties of 
ethanol. Molecular simulation studies using empirical 
force fields have addressed the equation of state, ther- 
modynamics, and structure and dynamics of solvation 
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] of aqueous ethanol solutions. The 
following general picture of the aqueous solvation in di- 
lute solutions has emerged: the hydroxyl group partic- 
ipates in the hydrogen bonded network, while the hy- 
drophobic alkyl group is accommodated in the hydrogen 
bonded network of water molecules. The nature of hy- 
dration structure around the hydrophobic part of alco- 
hols is still a controversial subject. It has been suggested 
0, |, i i 0, 111 that hydrophobic solutes enforces 
the network of hydrogen bonded water molecules around 
it and decreases their mobility, a notion referred to as 
hydrophobic hydration However, recent exper- 

imental and computational studies have shown that for 
small alcohols the structure and dynamics of the water 
molecules in the solvation shell is almost identical to that 
of bulk water ||, ||, ||. 

Much less is known for the solvation of ethylene in 
water under ambient conditions. Experimental [ p6[ and 
theoretical js^, ^ work addressed the clathrate hydrates 
of ethylene in water. The isolated ethylene-water com- 
plex has been a subject of various studies. In the low- 
est energy configuration, the ethylene molecule forms a 
weak bond with water. An early ab initio study of Del 
Bene in 1974 |^ has characterized this interaction as a 
TT-H hydrogen bond of a water proton with the tt elec- 
trons of the C=C bond. Several experimental techniques 
have been applied to measure the strength of this in- 
teraction, such as the matrix isolation study of Engdahl 
and Nelander Q, the microwave spectra study of 
Andrews and Kuczkowski and the molecular-beam 
measurements of Peterson and Klemperer 1 43 . Recently 
Tarakeshwar et al. H, nB, nQ and Dupre and Yappert 
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have performed calculations on the ethylene-water 
complexes with high level ab initio methods. The inter- 
action is weak compared to a hydrogen bonds such as 
in the water dimer. The role of the tt-H bond in aque- 
ous solvation under ambient conditions is still an open 
question. 

Molecular simulation provides an approach to study 
the microscopic behavior of liquids complementary to ex- 
perimental studies. All molecular simulations studies of 
aqueous ethanol solutions up to now are based on em- 
pirical force fields that are designed to reproduce a selec- 
tion of experimental data. Obviously, molecular simula- 
tions based on these potentials do not provide a picture 
completely independent from experiment. Moreover, the 
reliability of the results at conditions that are signifi- 
cantly different from those where the potential was de- 
signed for, may be questionable. Density functional the- 
ory (DFT) based molecular dynamics (MD) simulation, 
such as the Car-Parrinello molecular dynamics method 
p8| , where the interactions are calculated by accurate 
electronic structure calculations, provides a route to over- 
come these limitations. This has been demonstrated in 
studies of liquid water ^ |5^ and aqueous solva- 
tion lH, H ll. Important advantages of DFT-MD over 
force- field MD are that it intrinsically incorporates polar- 
ization, that it accounts for the intra-molecular motion 
and therefore allows for a direct comparison with spec- 
troscopy of intra-molecular vibrations, and that it yields 
detailed information on the electronic properties, such as 
the energy levels of electronic states and the charge distri- 
bution. In a broader chemical perspective it is important 
to note that DFT-MD is capable to study chemical re- 
actions in solution, where force-field MD would fail com- 
pletely as it cannot account for the change in chemical 
bonding. 

Here, we report on a DFT-based MD simulation of the 
solvation of ethanol and ethylene. First we describe the 
simulation methods. Then we show results of geometries 
and energetics of relevant gas-phase complexes, that will 
serve as a validation of the numerical methods employed. 
Subsequently, results of structure, dynamics, and polar- 
ization of the solvated species will be shown. We conclude 
with a discussion. 



II. COMPUTATIONAL METHODS 

Electronic structure calculations are performed using 
the Kohn-Sham H formulation of DFT H]. We em- 
ployed the gradient-corrected BLYP fimctional p8| , 
that has proven to give a good description of the struc- 
ture and dynamics of liquid water |Q . The DFT-based 
MD simulations of aqueous ethanol and ethylene are per- 
formed using the Car-Parrinello method as implemented 
in the CPMD package Norm-conserving Martins- 

Trouillier pseudopotentials are used to restrict the 
number of electronic states to those of the valence elec- 
trons. Cut-off radii for H, O and C atoms were chosen 



to be 0.50, 1.11 and 1.23 a.u., respectively both for the 
l=s and l=p terms. The Kohn-Sham orbitals are ex- 
panded in a plane wave basis, matching the periodicity 
of the periodic box with waves up to a kinetic energy of 
70 Ry. With this basis set energies and geometries are 
converged within 0.25 kcal/mol and 0.01 A , respectively. 
Vibrational frequencies are converged within 1 %, except 
for C-0 and 0-H stretch modes that are underestimated 
by 3 % and 5 % compared to the basis set limit values 
p3[ . In the molecular dynamics calculations, the ficti- 
tious mass associated with the plane-wave coefficients is 
chosen at 900 a.u., which allowed for a time step in the 
numerical integration of the cquations-of- motion of 0.145 
fs. 

To validate the computational approach we compared 
CPMD results for energies and structures of relevant 
gas-phase molecules and complexes with state-of-the-art 
atomic-orbital DFT-BLYP calculations performed with 
ADF |6^, ^ , and other high-level quantum chem- 

ical results taken from literature. The gas-phase calcula- 
tions with CPMD were performed in a large periodic box 
of 10 A using the screening technique of ||65| to elimi- 
nate the interactions among periodic images. We have 
not included zero-point energies in the energies of the 
gas-phase compounds. This also holds for computed en- 
ergies taken from literature and referred to in the present 
paper. This, to ensure a proper comparison between our 
results and those from literature. 

The MD simulations of the solutions were performed 
for a 'small' and a 'large' system to asses the finite size 
effects. For the small system a simulation of 10 ps was 
performed in a cubic periodic box of length 10.07 A, both 
for an ethanol solution of 31 waters and 1 ethanol as for 
an ethylene solution of 32 water and 1 ethylene. For the 
large system 5 ps simulations were performed using a pe- 
riodic box with bcc symmetry and a volume of 1977.6 
A3. This periodic cell, a truncated octahedron, is in 
shape closer to a sphere than a simple cube, and there- 
fore better suited for liquid simulation. The large ethanol 
solution consisted of 63 waters and 1 ethanol, while the 
large ethylene solution consisted of 64 waters and 1 ethy- 
lene. The box sizes for both the ethanol and ethylene 
solutions were set to match the experimental densities 
of the ethanol solutions under ambient conditions. For 
the ethylene system this will be slightly larger than the 
experimental density, as the effective volume of ethanol 
is a bit smaller than the combined volume of ethylene 
+ water. However, we do not expect this to give rise 
to observable changes in the calculated properties. For 
reference we performed 10 ps MD simulations of a sin- 
gle ethanol and ethylene in a periodic cubic box of 10 
A and a pure water systems of 32 water molecules in a 
cubic box of 9.85 A for simulation times of 10 ps. For all 
simulations there was an initial equilibration trajectory 
of 1 ps. Temperature was controlled by a Nose-Hoover 
thermostat ^ and fixed at 300 K. 
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III. GAS-PHASE COMPLEXES 

The ethanol monomer has two stable conformers very 
close in energy, the symmetric trans structure and the 
a-symmetric gauche structure (see fig. |^) . The main dis- 
tinction is the orientation of the OH bond with respect 
to the ceo plane. A microwave study has shown 
that the trans form is slightly (0.12 kcal/mol) more sta- 
ble than the gauche form. Table || lists the most impor- 
tant geometric data of the trans and gauche conform- 
ers, and compares the CPMD results with ADF calcu- 
lations, B3LYP calculations |70 , and experimental data 
0. CPMD and ADF bond lengths differ at most 0.01 
A and angles are within 0.5°. Comparing CPMD with 
B3LYP and experimental values yields differences upto 
0.03 A and 0.5°, and 0.04 A and 1° respectively. The 
calculated energy difference between the two conformers 
is listed in fig. |. CPMD, ADF, B3LYP of Ref. (7|, and 
MP4 |73) predict the trans conformer to be stable by 
0.07-0.10 kcal/mol, in good agreement with the experi- 
mental value of 0.12 kcal/mol. In contrast, the B3LYP 
calculation of ||7^ yields the opposite, with the gauche 
conformer stable by 0.23 kcal/mol. Vibrational frequen- 
cies obtained by a Fourier transform of the velocity auto 
correlation function (VACF) of a single ethanol at 300 K 
yields an OH stretch frequency of 3200 cm^^ and a CH 
stretch frequencies in the range of 2600-2800 cm~^. This 
should be compared to the experimental values of 3653 



cm-^ 721, and 2800-3000 em- 
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respectively. 



The tendency of BLYP to underestimate frequencies by 
« 10% is a known feature, and also observed in BLYP 
calculations of water Q and methanol Q . 

For the ethanol-water complex we distinguished four 
complexes, with both the ethanol trans and gauche con- 
former acting as proton donor or acceptor in the hydro- 
gen bond with water. Table Il| lists the most important 
geometric data obtained with CPMD and ADF. Fig. ^ 
shows the energy differences for all complexes. As for 
the ethanol monomer there is excellent agreement be- 
tween CPMD and ADF results, with deviations within 
0.02 A and 2° for bond lengths and angles, and the bind- 
ing energies within 0.2 kcal/mol. This indicates a state- 
of-the art accuracy of the electronic-structure method 
employed in CPMD. The CPMD-BLYP calculation pre- 
dicts the complex with the water molecule donating a 
proton to the ethanol gauche conformer to be the most 
stable, with a binding energy relative to the isolated wa- 
ter and trans ethanol conformer of 4.75 kcal/mol. The 
complex with the ethanol trans conformer donating a 
proton is less stable by 0.82 kcal/mol. Switching to the 
other ethanol conformer within either complex destabi- 
lizes the complex by 0.07 and 0.4 kcal/mol, respectively. 
The relative stability of the ethanol donor and acceptor 
complexes is similar to that for the methanol- water com- 
plex where CPMD-BLYP and the complete-basis-set 
MP2 estimate |7^ yields the methanol acceptor complex 
stable by 0.74 and 0.35 kcal/mol, respectively. There are 
no experimental data for the structure and energetics of 



the ethanol-water dimer. The MP2 result of Ref. |7|1 
provides the only high-level quantum mechanical study 
reported in literature. This study does not distinguish 
between the trans and gauche conformers. Therefore, a 
comparison with our results is somewhat limited. For 
the ethanol acceptor configuration the MP2 and CPMD- 
BLYP result for the complexation energy are similar. 
However, for the ethanol donor configuration the MP2 
complexation energy is more than 1.3 kcal/mol larger 
than the CPMD-BLYP results. Consequently, the MP2 
calculations yield an opposite relative stability of the two 
water-ethanol configurations, with an energy difference of 
0.92 kcal/mol. It should be noted that in the same study 
the methanol- water complex with the methanol as hydro- 
gen bond donor is found to be the most stable, in contra- 
diction with the MP2 basis-set limit result of Ref. [ [74[ . 
This suggests that in the Ref. a limited basis set 
or other factors could have lead to spurious reversal of 
the relative stability of the two water-ethanol configu- 
rations. The CPMD-BLYP hydrogen bond interaction 
energy for the ethanol-water dimer is of the same order 
of magnitude as the CPMD-BLYP result for the water- 
water [^^, water-methanol psf , and methanol-methanol 
[ |76| dimer. Comparison of these three dimers against 
high-level quantum chemical calculations and experimen- 
tal values indicates that BLYP underestimates the bind- 
ing energy by approximately 1 kcal/mol. For the ethanol- 
water dimer we could expect a similar difference. Here 
we should add that going from water, via methanol, to 
ethanol the dispersion forces become increasingly impor- 
tant. These are not accounted for in gradient-corrected 
functionals such as BLYP. Correlated methods such MP2 
incorporate, to a good approximation, dispersion forces. 
In Ref. ||7^ we estimated, in a comparison of BLYP and 
high-level MP2 calculations |^ for the methanol-dimer 
in non-hydrogen bonded configurations, that the absence 
of dispersion interaction in BLYP amounted to an under- 
estimate of the binding energy of wl kcal/mol. For the 
ethanol-water dimer this number could serve as an un- 
derestimate. Yet, although by far no insignificant, the 
magnitude of the deviation is much smaller than the 
hydrogen-bond interaction. Therefore, it can be argued 
that for a study of aqueous ethanol neglecting the dis- 
persion interaction is acceptable. 

Next we discuss the ethylene- water complexes. We 
will consider complexes with a single water and a wa- 
ter dimer. The dominant interaction of water with ethy- 
lene is a TT-H bond, where a proton of the water molecule 
binds to tt electrons of the double C=C bond. There ex- 
ist two stable complexes between ethylene and a single 
water molecule, indicated as the EWla and EWlb ge- 
ometry (see fig. Both geometries have Cs symmetry 
with the plane of the water molecule orthogonal to that 
of ethylene molecule. In the EWla structure the water- 
plane is parallel to the C=C bond, whereas in EWlb it 
is orthogonal. In table [II we compare the CPMD-BLYP 
geometries with ADF-BLYP and MP2[||, |^ and CCSD 
p6[ results. Fig. shows the energy differences, and 
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compares these to vahies reported in hterature. Again 
the CPMD and ADF results are in excellent agreement, 
with energies smaller than 0.02 kcal/mol, and geome- 
tries within 0.01 A and 0.3°, except for the non bonded 
OH„f, distance that is related to water orientation, a co- 
ordinate along with the energy surface is relatively flat. 
The BLYP result shows a significant tt-H binding en- 
ergy of « 1.4 kcal/mol. However, in MP2 j44[ |4|, Q 
and CCSD Q calculations, binding energies are ~ 1- 
2 kcal/mol stronger, accompanied by a shorter tt-H bond 
with differences of 0.12 A and 0.14 A for EWla and 
EWlb, respectively. The comparison with the experi- 
mental value gives a similar picture with the BLYP un- 
derestimating the binding energy « 2 kcal/mol. 



Matrix isolation studies ^ have revealed that the 
complex of ethylene with two water molecules consists of 
a water dimer that has one water molecule tt-H bonded 
to the C=C double bond. Recent high-level MP2 calcula- 
tions have indicated that the presence of the second water 
molecule enhances the strength of the tt-H bond^, ^t\ . 
From fig. ^ we see that BLYP result arc qualitatively 
in line with this observation, with a binding energy of 
the (water-dimer)-ethylene complex of 2.25 kcal/mol, up 
from 1.41 kcal/mol for the single-water binding. How- 
ever, quantitatively MP2 and BLYP compare less well, 
with BLYP underestimating both the total binding en- 
ergy as well as the increase from the single-water bind- 
ing. Again, the experimental value|^, |4^, that is closer 
to the MP2 result, shows a significantly stronger binding 
than BLYP. Both for the water-ethylene and the (water- 
dimer)-ethylene complex the differences of BLYP with 
MP2 and experiment are significant relative to the total 
binding energy. This absence of dispersion interaction in 
BLYP will be an important factor contribution to this 
discrepancy. 



Overall, we conclude that the reference calculations 
provide confidence that BLYP is capable of a quantita- 
tive study of a the aqueous solvation of a single ethanol 
molecule where interactions are dominated by relatively 
strong hydrogen bonds. BLYP qualitatively accounts 
for the weaker tt-H binding in water-ethylene systems. 
However, comparison with MP2 and experimental data 
suggest that its strength is significantly underestimated. 
Yet, for a single solvated ethylene the hydrogen bonds 
among the solvating water molecules will be the domi- 
nant interaction. Hence, we believe that BLYP will be 
able to quantitative described dilute aqueous ethylene 
solution. Both for ethanol and ethylene the absence of 
dispersion attraction in BLYP will have some impact, in 
particular for the coordination around the CH2 and CH3 
groups, where BLYP only accounts for the steric repul- 



IV. SOLVATION STRUCTURE 

A. Ethanol Solvation 

Fig. ^ shows ethanol-water radial distribution func- 
tions (RDFs) of the small and large ethanol solutions. 
The pronounced structure in the hydrogen bonding RDFs 
(HH, OH, HO, 00) are a clear indication of the presence 
of hydrogen bonds. All RDFs show that the small sys- 
tem gives a good description of the first solvation shell, 
while the large system also includes the second solva- 
tion shell. The peak positions of the force-field simu- 
lations of Fidler and Rodger indicated by crosses, 
are close to our results. Note that the good agree- 
ment for the CO RDFs, that are potentially sensitive 
to a proper description of the dispersion attraction, sug- 
gests that the absence of the dispersion interaction in our 
CPMD-BLYP simulation is of limited importance for a 
proper description of the aqueous solvation of the hy- 
drophobic group of ethanol. The position of the first 
peak of the OH- and HO-RDF indicates that the average 
hydrogen-bond length is 1.7 A for both the ethanol- 
donor and -acceptor bond. Integration of the 00-RDF 
upto the first minimum r — 3.3 A yields on average 3 wa- 
ter molecules in the first solvation shell of the hydroxyl 
group, in good agreement with the neutron diffraction 
value of « 3 ||2j] . Integrating the methyl-oxygen (C2O) 
RDF up to r = 5.7 A indicates that the first solvation 
shell of the methyl group consists of approximately 21 
water molecules. These coordination numbers are of the 
same order as the experimental estimation of Petrillo et 
al. who found that there are 18 ± 2 water molecules 
within 4 A from the center-of-mass of an ethanol in aque- 
ous solution. 



B. Ethylene solvation 

In fig. ^ we show the ethylene- water RDFs of the small 
and large ethylene solution. Both the carbon-oxygen and 
the intermolecular carbon-hydrogen RDFs are shown. 
Experimental or simulation data for ethylene in aqueous 
solution at room temperature were not available. The 
comparison between the small and large system shows a 
larger deviation in the first solvation shell than was found 
for the ethanol solutions. Apparently, the small ethylene 
system is not able to accommodate properly the first sol- 
vation shell. Integration the CO-RDF of the large system 
up to the first minimum at 5.7 A yields a hydration shell 
of 23 water molecules. More spatial information can be 
extracted from fig. ^ where, for the large system, the 
distribution of water H-atoms around the ethylene C=C 
axis is shown. X and Y are the components of the vec- 
tor joining a water hydrogen and the midpoint of the 
C=C bond. Here, X is the distance orthogonal to the 
C=C axis and Y the parallel distance. Note that we 
used the four-fold symmetry to improve the statistical 
accuracy yielding four identical quadrants. The ethylene 
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hydrogens are shown to visuaUze the size of the ethy- 
lene molecule, but do not indicate any angle dependence 
around the C=C axis. The figure shows a well defined 
elliptic region with no hydrogens present except for two 
weak, but clearly visible, peaks in the mid plane Y — 
a.t X ^ ±2.5. These should be attributed to the presence 
of vr — _ff bonded configurations. 



C. Water Structure 

In fig. the calculated water oxygen-oxygen RDFs of 
both the ethanol and the ethylene solutions are compared 
to those calculated for pure water. Note that the first 
solvation shell contains a large fraction of the total num- 
ber of water molecules, especially for the small system. 
Structural changes due to the solute molecule should 
therefore be detectable by comparing these systems with 
the pure water system. For the ethanol solution we see 
a small drop of the first peak accompanied by a slight 
increase at the first minimum. This indicates some de- 
crease in the hydrogen-bond structure, when compared 
to the pure water. As the small ethylene-water solu- 
tion is too small to accommodate a fully relaxed water 
solvation shell fig. Q only shows the RDFs of the large 
ethylene system. The RDFs show a similar behavior as 
for the ethanol solvation, with a slight decrease of the 
first peak and a small increase of the RDF in the re- 
gion around the first minimum. The small changes in 
the structure of the solvating water shell around ethanol 
and ethylene indicate that the hydrogen-bonded water 
network is very flexible and can easily open up to accom- 
modate small apolar solute groups without significantly 
changing its structure. Our findings are consistent with 
the neutron diffraction experiments of Turner and Soper 
p3| and force-field molecular dynamics simulation of Fi- 
dler and Rodger |^ who also did not find any evidence of 
structural enhancement in the hydration shell of ethanol. 
Turner and Soper found that only for larger hydropho- 
bic solutes this effect was experimentally detectable, but 
even then very small. The same trend was found for in a 
CPMD-BLYP molecular dynamics study of the solvation 
of methanol [^ . 

D. Hydrogen Bonds 

To examine the hydrogen bond statistics in the ethanol 
solution we adopted the definition of Ref. two 
molecules are hydrogen bonded if simultaneously the 
inter-oxygen distance is less than 3.4 A and the OHO 
angle is smaller than 30°. With this definition wc found 
that, for the large system, ethanol oxygen donates on 
average 0.9 hydrogen bonds and accepts 1.7. For the 
small system we found 1.0 and 1.6, respectively. This in 
consistent with the fact that approximately three water 
molecules occupy the first solvation shell of the hydroxyl 
group. For comparison, in a CPMD-BLYP simulation 



we found that methanol in a dilute aqueous solution do- 
nates on average 0.9 hydrogen bonds and accepts 1.5 
pSf . ^From the pure water simulation these numbers 
were measured to be 1.7. 

For the ethylene molecule a well defined definition for 
the TT-H hydrogen bond does not exist. To investigate 
the infiuence of the tt-H hydrogen bond on the solva- 
tion structure we looked at the water hydrogen positions 
relative to the C=C axis (see fig^). Inside the elliptic re- 
gion around the C=C bond that is depleted of hydrogens 
we clearly detect near the center-of-mass of the ethylene 
molecule at a distance of « 2.5 A a region with an in- 
creased amount of water hydrogens. Integration over this 
region with < X < 2.6 and < y < 1 yields a value of 
0.42. This implies that approximately 40% of the time 
a water molecule is oriented towards the double bond, 
forming a tt-H bond. 



V. SOLUTE DYNAMICS 

The time scale of the present simulations (5-10 ps) 
allows for an analysis of the short-time dynamics. As 
we mentioned above BLYP underestimates most vibra- 
tional frequencies by more than 10 %. However, the fre- 
quency shifts upon solvation still provide valuable infor- 
mation and allow for a direct comparison to experiments. 
Fig. ^ shows the spectrum of the velocity auto correla- 
tion functions (VACFs) of the ethanol hydroxyl hydrogen 
for the large ethanol solution and for the single isolated 
ethanol. The spectrum is of limited accuracy due to the 
relative short trajectory. However, the calculations show 
that, upon solvation, there is a clear red shift of about 
200 cm^^ of the OH stretch frequency. Experimental 
data for the OH frequency shift of ethanol in dilute aque- 
ous solution is not available. However, the OH red shift 
is a typical characteristic of a hydrogen bonded liquid. A 
comparison with measured shifts in liquid ethanol, from 
3676 cm~^ ||7^ to 3330 cm^^Q, shows a similar trend. 
In a CPMD-BLYP simulation of a dilute aqueous solution 
of methanol we found a similar red shift for the methanol 
OH frequency of about 200 cm^^ 

The vibrational spectrum of the VACF of the ethylene 
molecule is not shown. As for the ethanol system there 
was a limited accuracy due to the relatively short calcu- 
lated trajectory. We observed that upon solvation the CH 
and CC peaks do broaden. However an estimate for peak 
shift falls outside the accuracy of the calculated spectra. 
The matrix isolation studies of Engdahl and Nelander 
pO| , ^ show minor changes in frequencies when isolated 
ethylene is compared with the ethylene-water complex, 
with the largest shift being a blue shift of 12 cm^^ of the 
out-of-plane bending mode {vr w 947 cm"-'^). 
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VI. POLARIZATION 

As the electronic structure is an intrinsic part of a 
CPMD simulation, detailed information about the elec- 
tronic charge distribution can be obtained. To quan- 
tify the charge distribution we used the method of max- 
imally localized Wannier functions that transforms the 
Kohn-Sham orbitals into Wannier functions, whose cen- 
ters (WFC) can be assigned a chemical meaning such as 
being associated with an electron bonding- or lone-pair 
1^ , p2| . It also provides a unambiguous route to de- 
termine the dipole moment of individual molecules in a 
condensed phase by assuming the electronic charge to be 
distributed as point charges located on the WFCs. 

Table |^ list the calculated dipole moments for the 
gas-phase systems and the average dipole moment for 
the solvated systems. The latter were obtained from 28 
independent configurations of the small solvated systems. 
The distributions of these dipole moments are shown in 
fig. ^. For ethanol we observe an significant increase in 
going from the gas-phase molecule via the water-ethanol 
cluster to the fully solvated system. This trend is very 
similar to what is found is found for water [ p3| , p4[ and 
methanol fz^ and may be considered typical feature of 
a strongly hydrogen bonding molecule. The calculated 
distribution of dipole moments of solvated ethanol shows 
that thermally driven fluctuations give rise to a signifi- 
cant variation ranging from 2.0 D to 4.0 D. 

Also for ethylene we observed a significant change of 
the dipole moment upon solvation. Being apolar in the 
gas phase, the average dipole moment increases up to 
0.5 D when complexed with the water dimer or when in 
solution. In solution it exhibits rare fluctuations where 
the dipole moment reaches values of up to 1.0 D. Fig. |l^ 
shows a typical snapshot where such an extreme high 
value of the dipole moment is reached. The two WFCs 
of the ethylene double bond are located near the mid- 
dle of the C=C bond, just below and above the plane 
of the molecule. In the figure we see clearly that the 
upper TT-WFC is acting as proton acceptor. The tt-H 
bond shifts this WFC further out of plane inducing an 
ethylene dipole moment orthogonal to the plane of the 
molecule. The snapshot suggests an enhancement of the 
ethylene polarization by the fact that the tt-H bonded wa- 
ter molecule has in turn a rare four-fold hydrogen-bonded 
water coordination, with three of the waters donating a 
proton. This should be seen as a manifestation of the 
strengthening of the tt-H bond upon complcxing water 
molecules with the water molecule that is tt-H bonded to 
the ethylene, a feature found in gas-phase MP2 calcula- 
tions ^ and also in the present BLYP calculations, 
and discussed above in section III. 



VII. CONCLUSIONS 

We have studied the solvation of ethanol and ethylene 
in water by DFT-based (Car-Parrinello) Molecular Dy- 



namics. Validation of the computational approach by 
comparing structure and energetics of relevant gas-phase 
complexes against experimental results and state of the 
art quantum chemical calculations showed that CPMD 
employing the DFT-BLYP is capable of qualitatively de- 
scribing the aqueous solvation of a single ethanol of ethy- 
lene molecule. 

The structural properties of the ethanol solvation were 
in good agreement with both neutron diffraction data 
p3[ and force-field MD simulations |Q . We found that in 
aqueous solution ethanol accepts on average 1.7 hydrogen 
bonds and donates on average 0.9 hydrogen bonds. For 
ethylene we found it has approximately 0.4 tt-H bonds 
with a water molecule. Both for ethanol and ethylene 
the simulations provided no structural evidence for hy- 
drophobic hydration: the structure of the hydrating wa- 
ter shell was not enhanced compared to that of pure wa- 
ter. The calculation even indicated a slight decrease in 
the structure. For aqueous ethanol the calculated red 
shift of the hydroxyl vibration upon solvation was con- 
sistent with experimental findings. 

Analysis of the electronic charge distribution by means 
of Wannier functions showed that the interaction with 
water can significantly increase the dipole moment of 
the ethanol and the ethylene molecule. The average 
dipole moment of ethanol increases from 1.8 D in the 
gas phase to 3.1 D in aqueous solution. Ethylene, that 
is apolar in the gas phase, attains an average dipole mo- 
ment of 0.5 D in solution. We identified configurations, 
with a TT-H bonded water molecule that has a rare four- 
fold hydrogen-bonded water coordination, where the in- 
stantaneous dipole moment of ethylene takes values of 
up to 1 D. Such configurations with large solute dipole 
moments may also play an important role in activating 
chemical reactions involving the solutes as we have seen 
in a CPMD-BLYP study of the acid-catalyzed hydra- 
tion of ethene |jl^. The strong polarization effect raises 
questions towards the common consideration of think- 
ing ethylene as a apolar and hydrophobic molecule. The 
electronic charge analysis also points out the necessity 
of polarizable force fields for both ethanol and ethylene 
when dissolved in water. The present results may be 
considered valuable for the design of such force-fields. 
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TABLES 

TABLE I: Gas-phase complexes; ethanol gauche and trans monomer. Distances in A, angles in degrees. Our results: CPMD- 
BLYP and ADF-BLYP, are compared with B3LYP/6-311G(2d,2p) and experimental microwave spectroscopy data fri\. 



Geometry 


Gauche 


Method 


CPMD 


ADF 


B3LYP 


EXP 


^■OH 


0.981 


0.971 


0.961 


0.945 


'"CO 


1.455 


1.447 


1.429 


1.427 


'"CC 


1.529 


1.5S2 


1.521 


1.530 


< rcH > (in CHS) 


1.099 


1.097 


1.091 


1.094 


< rcH > (in CH2) 


1.100 


1.099 


1.092 


1.094 


QCOH 


108.0 


108.2 


108.7 


108.3 


QCCO 


IIS.O 


113.1 


113.0 


112.2 


Geometry 


Trans 


Method 


CPMD 


ADF 


B3LYP 


EXP 


roH 


0.980 


0.970 


0.960 


0.945 


'"CO 


1.458 


1.449 


1.432 


1.425 


'"CC 


1.52S 


1.524 


1.515 


1.530 


< rcH > (in CHS) 


1.098 


1.096 


1.090 


1.094 


< rcH > (in CH2) 


I.IOS 


1.102 


1.095 


1.094 


QCOH 


108.4 


108.4 


109.0 


108.3 


QCCO 


107.6 


108.0 


108.0 


107.2 



TABLE H: Gas-Phase complexes: ethanol-water dimers. Distances in A, angles in degrees. CPMD-BLYP and ADF-BLYP are 
compared. We differentiate four complexes with the ethanol can be in the gauche or trans geometry and acts as proton acceptor 
or proton donor in the ethanol-water hydrogen bond. H„6 is the hydrogen of the water that is not involved in a hydrogen bond, 
Ci is the hydroxyl carbon, C2 is the methyl carbon and O^, is the water oxygen. 



Geometry 


Gauche 




Trans 




H-bond 


H-acc 


H-don 


H-acc 


H-don 


Method 


CPMD 


ADF 


CPMD 


ADF 


CPMD 


ADF 


CPMD 


ADF 


roH in H20 


0.994 


0.981 


0.981 


0.972 


0.993 


0.982 


0.981 


0.972 


roH„, in H20 


0.981 


0.970 


0.981 


0.971 


0.978 


0.969 


0.982 


0.971 


roH in eth. 


0.981 


0.971 


0.987 


0.978 


0.982 


0.970 


0.988 


0.977 


roH H-bond 


1.934 


1.912 


1.988 


2.011 


1.928 


1.929 


1.960 


1.993 


'"CC 


1.526 


1.528 


1.532 


1.533 


1.521 


1.522 


1.523 


1.526 


rco 


1.469 


1.457 


1.452 


1.440 


1.474 


1.460 


1.454 


1.442 


roo 


2.924 


2.884 


2.978 


2.983 


2.908 


2.989 


2.988 


2.968 




3.653 


3.626 


3.773 


3.793 


3.641 


3.641 


3.815 


3.764 


ro„C2 


3.932 


3.920 


4.128 


4.033 


4.013 


4.008 


5.229 


5.189 


a HOH in H20 


104.2 


105.0 


104.6 


104.9 


104.6 


105.0 


104.2 


104.9 


a COH 


107.9 


108.4 


108.5 


108.7 


107.9 


108.2 


108.9 


108.5 


Q OH-O 


172.6 


170.2 


172.4 


172.8 


168.1 


168.4 


175.0 


174.1 
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TABLE III: Gas-Phase complexes: ethylene- water, EWla, EWlb structures (see fig. Distances in A, angles in degrees. 
CPMD-BLYP, ADF-BLYP and MP2-TZ2P-I-I- Q calculations are compared. H is the water-hydrogen involved in the vr-H 
bond; H„i, is the other hydrogen of the water molecule. CM is the midpoint of the two carbons of the ethylene molecule. 



Geometry 


EWla 


EWlb 


Method 


CPMD 


ADF 


MP2 


CPMD 


ADF 


MP2 




1.334 


1.336 




1.335 


1.336 






2.483 


2.485 


2.363 


2.570 


2.524 


2.383 


rO-CM 


3.452 


3.460 


3.301 


3.527 


3.497 


3.337 


?'H„6-CM 


3.909 


3.833 


3.816 


3.957 


3.881 


3.793 




0.984 


0.975 


0.962 


0.978 


0.974 


0.962 


''OH„i, 


0.980 


0.970 


0.958 


0.976 


0.970 


0.958 


qhoh 


104.2 


104.5 


104.7 


104.4 


104.5 


104.7 



TABLE IV: Dipole moments obtained by Wannier-function analysis. The liquid water value (a) was taken from Silvestrelli et 
al. |84| . In the last column the total dipole moments of the gas-phase complexes are given. The first value for water in 
the E2W complex is the water closest to the ethylene, the second value is the dipole of the other water molecule. The solvated 
ethanol and ethylene dipole moments are obtained by taking the average of 28 independent configurations from the calculated 
trajectories. 



complex 


ethanol / ethylene 


water 


total 


single- water 




1.82 


1.82 


liquid water 




3.0 (a) 




ethanol gauche 


1.83 




1.83 


ethanol trans 


1.66 




1.66 


gauche P-acc 


2.32 


2.16 


2.76 


gauche P-don 


1.91 


1.91 


2.62 


trans P-acc 


2.13 


2.13 


2.47 


trans P-don 


1.91 


2.12 


2.29 


solvated ethanol 


3.08 


2.97 




ethylene 


0.00 




0.00 


EWla 


0.37 


1.95 


2.10 


EWlb 


0.33 


1.90 


2.02 


E2W 


0.51 


2.25 , 2.17 


2.41 


solvated ethylene 


0.51 


3.00 
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FIGURES 






FIG. 1: Illustration of the trans and gauche conformers of ethanol. 



^ ethene ^ ^ 



1-1 water 



t S™''"^ )2 

trans ^ 

ethanol ^ 



+1 water \ 



EWlb 
EWla 



, + 1 water 



E2W ' 



gauche H-don 

i 

trans H-don ^ i 

trans H-acc ^ 5 

gauche H-acc 



This work 


literature 


1 


CPMD 


ADF 


Computational 


Exp. 


1 


9.09 


9.35 






2 


0.08 


0.08 


0.07 (a), 0.10(b), -0.23 (c) 


0.12(d) 


3 


4.75 


4.60 


4.4-5.4 (e) 




4 


0.07 


0.09 






5 


0.82 


0.81 


-0.92 (e) 




6 


1.22 


1.00 






7 


1.45 


1.41 


2.3-2.9 (f,g), 2.0-3.3 (h) 


3.7 (i) 


8 


0.05 


0.06 


0.13 (f) 




9 


2.27 


2.25 


3.8-5.6 (j), 3.8-4.7 (g) 


5.5-6.1 (k) 



FIG. 2: Relative energies of ethylene-water and ethanol-water in kcal/mol. Zero-point energies (ZPE) are not included. For 
ethylene, the energies are relative to the separate molecules, whereas for ethanol the separate water and trans ethanol conformer 
is taken as the reference value. 3 and 7 indicate the binding energies of ethanol and ethylene respectively. 9 is the binding energy 
between ethylene and the water-dimer. The literature values are: a: MP4-(SDQT)/cc-pVTZ b: DFT/B3LYP/6-311G(d,p) 
[g, c: DFT/B3LYP/6-311G(d^ d: microwave experimentj| e: MP2/6-311+G(2df,2pl+BSSE (75|, f: MP2/TZ2P++ 
ig, g: MP2/6-311+G(2d,2p) GtI, h: CCSD(T)/aug-cc-pVDZ ^f^i: Matrix isolation study Q, j: MP2/aug-cc-pVDZ 
and k: Matrix isolation study |40[ . For (e), (f), (g), and (h) we gave the BSSE-corrected (lowest value) and the non-BSSE 
corrected values (highest value)^i) and (k) are the experimental values minus the ZPE of |47[. 



f > 

FIG. 3: Illustration of the EWla and EWlb ethylene-water complexes. 




FIG. 4: Ethanol-water RDFs of the ethanol solutions. All insets show both the results of the small (solid line) and the large 
system (dashed line) . For all graphs gAB indicates that the first atom A belongs to ethanol and the second atom B belongs to 
water. If the first is a hydrogen, A=H, then always the hydroxyl hydrogen of the ethanol is meant. Ci is the carbon bonded 
to the hydroxyl group, C2 is the carbon of the methyl group. 
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1.6 
1.2 
0.8 
0.4 

0.0. 





gCH 



1.0 2.0 3.0 4.0 5.0 6.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 

r r 

FIG. 5: Ethylcne-water RDFs of the ethylene solutions. The result for small system and large system are given by the solid 
and dashed line, respectively. H denotes the water hydrogens only. 



8.0 



4.0 



o< 



0.0 



-4.0 



-8.0. 




-8.0 



-4.0 



0.0 



4.0 



8.0 



X(A) 



FIG. 6: Time and rotational averaged normalized density distribution p{r)/po (with po is the average hydrogen density) of 
the water-hydrogens around the solvated ethylene molecule obtained from the large system simulation. The rotational average 
is about the C=C axis. The four-fold symmetry is imposed so that all quadrants contain the same information. The density 
regions 1.0 < p{r)/po < 2.0 and p{r)/po > 2.0 are illustrated by small and large pixels, respectively. 
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FIG. 7: Water-water RDFs. Left: the two ethanol solution systems are compared with pure water. Right: the large ethylene 
solution system is compared with pure water. Thin solid lines indicate the pure water, thick solid lines the small system, and 
thick dashed lines the large system. 




FIG. 8: Spectrum of the VACF of the hydroxyl ethanol hydrogens in the gas-phase and in the large aqueous system. 
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dipole moment (D) 

FIG. 9: Dipole moment distribution of the solvated ethylene and the solvated ethanol molecule. 




FIG. 10: A snapshot of ethylene solvated in water. The figure shows a rare configuration where the ethylene molecule has a 
large dipole moment of « 1 D. For clarity, only the ethylene molecule and five water molecules are shown. Carbon and oxygens 
arc dark grey. Hydrogens are light grey. Besides the atoms also the WCFs are shown (middle grey). The water molecules have 
WFCs in both OH bonds and two aside the oxygen corresponding to the lone pairs. Ethylene has four WFCs along the CH 
bonds and two out-of-plane WFCs between the carbons corresponding to the tt bond. The dashed lines indicate the hydrogen 
bonds. One the protons of the central water molecules points towards the 'upper' tt bond WFC. This is the manifestation of 
the TT-H bond in this picture. 



